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Previously a compact formula for total reaction cross section for heavy-ion collisions as a function 
of energy was obtained by treating the angular momentum I as a continuous variable. The accuracy 
of the continuum approximation is assessed and corrections are evaluated. The accuracy of the 
compact equation can be improved by a simple modification, if a higher accuracy is required. Simple 
rules to determine the barrier heights and the penetration probability for the I partial wave from 
experimental data are presented, for the collision of identical or non-identical light nuclei. 
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I. INTRODUCTION 

Nuclear fusion in heavy-ion collisions is an important 
process in many phenomena. The sub-barrier fusion of 
light nuclei plays an important role in the evolution of 
massive stars, the dynamics of white dwarf Type la su- 
pernovae, and explosions on the surface of neutron stars 
The fusion of heavy nuclei just above the barrier is 
an important tool in the production of superheavy nu- 
clei 0. The fusion of identical light nuclei at high en- 
ergies reveals interesting effects in isolated high-angular- 
momentum states of the fused system Q . Nuclear fusion 
of light nuclei is utilized in accelerator-based heavy-ion 
inertial fusion for fusion energy production 



Previously, a simple analytical expression was obtained 
for the total reaction cross section ay in the collision of 
nuclei A\ and Ai as a function of the collision energy E 




where Eq is the barrier height for the 1=0 partial wave, 
T=hiu/2iT is the energy width parameter in the potential 
barrier penetration probability, u is the frequency of the 
inverted parabola representing the potential barrier, and 
Rq = r o(^4-i^ 3 + A^ 3 ) is the spatial separation between 
the two nuclei at the potential barrier. By assuming that 
the fusion process is equivalent to the strong-absorption 
of ingoing waves passing through the potential barrier, 
the total reaction cross section ([I) can be interpreted as 
the heavy-ion fusion cross section. 

Even though the fusion process involves complicated 
dynamics of channel coupling, dynamical distortions, po- 
larizations, deformations, particle transfers, proximity 
interactions, and barrier penetrations @, the simple ex- 
pression of Eq. ([T]) provides an efficient way to repre- 
sent experimental data in terms of important physical 
parameters, whose systcmatics give valuable insights into 
the dynamics of the process. The relationship between 
the barrier penetration model used in Ref. [5| and the 
coupled-channels calculations for heavy-ion fusion was 
discussed previously Q. Equation ^ has been success- 
fully applied to describe fusion cross sections in the col- 
lision of many projectile and target combinations @ . 



As the range of fusion energy in astrophysical processes 
and fusion-energy production extends to the region be- 
yond the sensitivity of present-day measurements, theo- 
retical extrapolations are needed to access many relevant 
rates of fusion reactions [l|. It is desirable to examine 
the accuracy of the simple formula ([1]) in the sub-barrier 
region. In the other extreme in higher energy fusion, the 
recent interesting observation of the stepwise increase of 
the fusion cross section for two identical spin-0 nuclei Q 
also calls for an analysis on its accuracy in the higher 
energy regime. 

The simple result of Eq. (1) was obtained by treating 
the discrete angular momentum I as a continuous vari- 
able. Such a treatment incurs errors. We therefore wish 
to assess the accuracy of Eq. ((T|) over the whole energy 
range and to find how it may be improved if a higher 
accuracy is desired. 

The simple result of Eq. (1) relies also on the assump- 
tion that the fusion barrier height Ei for the I partial wave 
is a linear function of 1(1 + 1). While such an assumption 
is a reasonable concept for many reactions, there are nu- 
clear collisions for which such an assumption is not valid, 
as is evidenced by the disagreement of the 12 C+ 13 C data 
[HHI with the main features of Eq. (1). To diagnose such 
a pathological case, it will be useful to develop tools that 
will enable us to determine the fusion barrier heights Ei 
and the penetration probabilities P\ (E) , as a function of 
I from experimental data. Furthermore, the direct de- 
termination of these physical quantities provides useful 
insight into the fusion process. 

For those cases where the basic assumption of the lin- 
ear dependence of the fusion barrier Ei on 1(1 + 1) does 
not hold, we would like to propose alternative methods 
to describe the reaction cross section. We wish to design 
a framework to partition the reaction cross section such 
that contributions to different regions of I partial waves 
can be singled out for special scrutiny. 

This paper is organized as follows. In Sec. II, we eval- 
uate the reaction cross section and develop the rules for 
fusion barrier analysis. In Sec. Ill, we present the for- 
mulation of the reaction cross section, the continuum ap- 
proximation, and its corrections. In Sec. IV, we give the 
numerical results and the comparison with experimen- 
tal data. In Sec. V, we carry out a barrier analysis for 
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12 C+ 13 C and show that the barrier E\ for that reaction 
is not a linear function of 1(1 + 1). We show how the 
fusion cross section of such a pathological case can be 
described by an alternative method of partitioning the 
reaction cross section. We present the rules for barrier 
analysis in Sec. VI and the rules for the penetration prob- 
ability analysis in Sec. VII, for the collision of identical 
or non-identical light nuclei. In Sec. VIII, we present our 
conclusions and discussions. 



II. REACTION CROSS SECTION AND FUSION 
BARRIER ANALYSIS 

We approximate various barriers for different par- 
tial I waves by inverted harmonic-oscillator potentials of 
height Ei and frequency uji in the ingoing- wave strong- 
absorption model. For a collision energy E, the probabil- 
ity for the absorption of the I partial wave is then given 
by the Hill- Wheeler penetration probability formula , 



Pi(E) 



1 



l + exp{2n(Ei - E)/tkdi}' 



(2) 



In consequence, the total reaction cross section as a func- 
tion of E for the collision of unequal nuclei is [f| 



21 + 1 



1=0,1,2... 



+ exp{2n(Ei - E)/hui}' 



(3) 



The above expression can be cast in a more illuminat- 
ing form in terms of the deBroglie wave length A = l/k. 
The natural unit of cross sectional area in deBroglie wave 
length scales is 7rA 2 [ll| , which can be conveniently called 
the deBroglie cross section. Using Eq. ([3]), we can con- 
struct the dimensionless measure of the reaction cross 
section oy in units of 7rA 2 , (cr r /7rA 2 ), at the collision en- 
ergy E, as given by 

(a r /ir% 2 )(E) = /(0 = fo + h + h + h + ••, (4) 



(=0,1.2, 



where 



fl = f(l) 



21 + 1 



i + cxp{(£, -£)/r z } ; 



(5) 



and Ti = hwi/2ir. We can evaluate (a r /-KX 2 )(E) at the 
collision energy E. We find from Eq. ^ that for the 
energy E such that 



I 2 < (a r /ir\ 2 )(E) < (Z + l) 2 , 
the reaction cross section at energy E is 

2Z + 1 



(a r /Trr)(E) = l z + 



1 + exp{(£, - E)/T) 



(6) 



C(l,E),(7) 



where the correction term C(l,E) takes into account the 
width Ti for the barrier penetration. It is given explicitly 



by 



C(l,E) 



i-i 



21' + 1 



l'=0 
oo 



i'=i+i 



exp{(E-E v )/T v } 
2V + 1 



e(z-i) 



exp{(Ei> -E)/r v y 



(8) 



where Q(x) = 1 for x > 0. 

For the evaluation of the reaction cross section and the 
correction term C(l, E), we shall study a simple model in 
which we assume that the barriers Ei and the frequencies 
fajJi (or Ti) are related to I by 



Ei = E + 



h 2 l{l + l) 



2fiR 2 ' 

hui ~ fiw, ( or r ; = r), 



(9) 
(10) 



where /i = AiA2m mlc \ eon / (A\ + A2) is the reduced mass. 
We shall further convert the summations in Eq. © as 
integrals in the continuum approximation, then the cor- 
rection term is given explicitly by 



C(l =0,E): 

C(l = 1,E) 



^ 2 ~ ln 



2 M i? 2 r ln 



n 2 



1 + exp{(£ - £d)/r} 

1 + exp{(E - E 2 )/T) 
1 



,(H) 



l + exp{(E-E )/ry 



(12) 



C(l > 2,E) 



2/ii^r ln 



ft 2 

2 M i? 2 r 
h 2 



in 



i + ex P {(£-£ i+1 )/r} 
i + cxp{(£; t _i -g)/r} 

l + exp{(£ -E)/T} 



.(13) 



Thus, Eq. (UJ), with supplementary equations ©, (fTT|) . 
(fT2|) . and (fl3]) . gives the reaction cross section as a func- 
tion of energy E. 

We can evaluate the reaction cross section 
(oy/7tA 2 )(£?/) at the barrier Ei. It is given by 



(a r /irX 2 ) 



= l(l + l) + ±+C(l,E l ), (14) 



where the correction term C(l,Ei) is 
2T 



C(l = 0,Ei) 
C(l = l,Et) 

C(l>2,Ei) 



Ei — Eq 

2T 

El — Eq 



l + exp{(El-Ei +1 )/T} 
l + exp{{Ei- E [+1 )/T} 



1 



l + exp{(Ei-Eo)/ry 
2T 



Ei — Eq 

2T 

Ei — Eq 



l + exp{{Ei-E l+1 )/T] 

l + exp{(£ t _i -Efi/T} 
l + expi(Eo-Ei)/T} 



(.15) 
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Equation (|14[) has a simple physical interpretation. As 
illustrated in Fig. 2.1 of Blatt and Wcisskopf [ll[, the 
partial wave I' contributes 21' + 1 units to the dimcn- 
sionless measure of the reaction cross section. The total 
contribution is the integral of J dl'(2l' + 1). Therefore, 
the dimensionless measure of the reaction cross section, 
up to the fusion barrier of the I partial wave, is given 
by 1(1 + 1) on the right hand side. The additional con- 
stant 1/2 is purely quantum mechanical in origin and it 
depends on the symmetry of the colliding system, as will 
be discussed in Sec. VI. The correction term C(l,Ei) 
in Eq. (fl"4| takes into account the finite energy width Ti 
for barrier penetration probability relative to the spacing 
between adjacent fusion barriers. 

Equations (fl"4|) and (fl~5|) can be inverted to provide 
the rule for the "barrier analysis" for unequal nuclei as 
follows. The fusion barrier Ei for the I partial wave is 
located at the value of energy E at which the dimen- 
sionless reaction cross section measure, oy/-7rA 2 , is equal 
to 1(1 + l)+l/2+C(Z, Ei). If the dimensionless measure 
(<7 r /7rA 2 ) can be obtained experimentally as a function 
of E, the heights of various fusion barriers Ei can be de- 
termined iteratively, within the present model of fusion 
barrier penetration. 

In the beginning of the iteration, one neglects the cor- 
rection C(l,Ei), and the barriers Ei can be determined 
from the ay/VA values. With the knowledge of the bar- 
riers heights Ei the correction terms C(l, E{) can be eval- 
uated for different I partial waves, and the barrier quan- 
tities Ei can be corrected. In these iterations, it is nec- 
essary to know the width parameter T, which can be ob- 
tained either from a fit of the experimental fusion cross 
section with Eq. ([T]) , or from the penetration probability 
analysis, as will be discussed in Sec. VII. 

We list the values of cr r /7rA 2 at which the fusion bar- 
riers Ei are located for the collision of unequal nuclei in 
Table I, when the correction term C(l,Ei) can be ne- 
glected. 



TABLE I. The value of the dimensionless measure a r /nX 2 at 
which the fusion barriers Ei for the I partial wave is located, 
for the collision of unequal nuclei, when the correction term 
C(l,Ei) can be neglected. 
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The ratio T/(Ei — Ei-\) in the correction term C(l,E{) 
varies with the colliding nuclei mass number as A b / 3 T /2l. 
Thus, a decrease in the mass number or an increase in the 
angular momentum I will lead to a smaller V / (Ei — -Ej-i) 
and a smaller correction term C(l,E{). Our investi- 
gations in subsequent sections (Tables II and III) in- 
dicate that for light nuclei collisions, the condition of 
r <C \Ei —Ei±i\ is approximately fulfilled so that C(l, E{) 
is small. It becomes appropriate to neglect the correction 
term in the barrier analysis for light nuclei collisions. For 
heavy nuclei collisions, while the neglect of the correction 



C(l,E{) may be appropriate in the barrier analysis for 
large angular momentum I, the correction term must be 
properly taken into account for partial waves with small 
values of I . 



III. THE CONTINUUM APPROXIMATION 
AND ITS CORRECTIONS 

Our model assumption of a linear dependence of Ei 
on 1(1 + 1) in Eq. ([9]) is a reasonable concept for cases 
when the effective separation of the two colliding nuclei 
Rq at the fusion barrier is insensitive to the change of 
the angular momentum I. While such an assumption is 
reasonable for most reactions, there are however cases, 
such as 12 C+ 13 C, in which such an assumption may not 
be valid. The method of barrier analysis we have just 
developed in Sec. II may be used to diagnose the patho- 
logical case. We shall discuss the collision of 12 C+ 13 C in 
Sec. V. 

By replacing the sum in Eq. ([3]) by an integral in the 
continuum approximation, the reaction cross section can 
be integrated to yield the analytical formula of Eq. ([1]) 
[U- Such a replacement of the discrete I variable by a 
continuous variable incurs errors. It is desirable to find 
the magnitude of the errors and ways to correct for these 
errors, if a higher accuracy is required. For brevity of 
notation, we introduce a and g to rewrite /; as 



fl = f(l) 



21 + 1 



21 + 1 



l + cxp{[aZ(Z + l)-e)]/r} 1+g 



,(16) 



where g = cxp{ [al(l + 1) — e] / T}, e = E — E , and 
a = h 2 /2fiRl 

Our effort to examine the errors brings us to parti- 
tion the contributions in Eq. (j4]) into two groups: (i) one 
group of I states for which the continuum approxima- 
tion is a reasonable concept and analytical results can 
be readily obtained, and (ii) another group of discrete 
I states which remain as they are, without applying the 
continuum approximation, and their contributions to the 
total reaction cross section can be subsequently singled 
out for scrutiny. 

The I = state is important in sub-barrier fusion and 
it is not suitable for the continuum approximation. We 
shall keep /o to remain as it is in Eq. (Q| . We can express 
fi with I > 1 as a continuous integral with a correction 
A/j 



W+l/2 

fi= / dlf(l) + Afi. 

Jl-X/2 

The function f(l) has an indefinite integral 



dlf(l)=— In ■ 
a 



exp 



al(l + 1) 



(17) 



=F(I)(18) 



In terms of the function F(l), we have 

/, =F(/ + l/2)-F(Z-l/2) + A/,. (19) 
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By definition, the correction term Afi is then given by 

A/, = f(l) - [F(l + 1/2) - F(l - 1/2)]. (20) 

Treating I as a continuous variable in the above equation 
and expanding the function F(Z±l/2) about I in a Taylor 
series with Al = 1/2, we obtain explicitly 



Afi = -2 



E 

n=2,4,. 



(AZ) n+1 d 



(n + iy. dl 



(21) 



We thus obtain the central result that in the continuum 
approximation, any term fi with I > 1 in the series of 
Eq. can be replaced by Eq. (|TT)|) of the known function 
-F(Z), with A/; given by Eq. ([2~Tj) . For example, if we wish 
to partition the partial wave into those with [0, li] as a 
discrete sum, with those in [1^ + 1, oo] in the continuum 
approximation, then we obtain for such a partition 



^ = J2f(l)-F(l L + 1/2) 




( n + iy.di nn > 



(22) 



We can use the relation 



is the correction due to the continuum approximation up 
to the third order in AZ=l/2. 

The above considerations can be generalized. For the 
most general case \l2\) . 



ttA 2 



£[i + ?/ (-i)'1/(U 



(27) 



where (i) J] = and l v = v = 0,1,2,3,... for unequal 
nuclei, (ii) rj = 1 and l v = 2v = 0,2,4,... for identical 
spin-0 nuclei or for identical spin- 1/2 nuclei with sym- 
metric spatial and antisymmetric spin wave functions, 
and (iii) rj = —1 and l v = 2v + 1 = 1, 3, 5, .. for identical 
spin-1/2 nuclei with antisymmetric spatial and symmet- 
ric spin wave functions. The sum over l v can be converted 
into a sum over v with v = 0, 1,2, ... We obtain up to the 
third order (Ai^) 3 with Av = 1/2, 



2 dl v 



2/n 



1 



dv 1 + cxp{[aZ (/ + 1) - e)]/r} 



+-kR z —\u\ 1 + exp 



h/2(l 



/2l'l/2 



l)a 



-2-kX 



dl v 
dv 



E 



^=1,2,3, 



r 

(A^) 3 d 2 
3! d& 



f(Q • (28) 



TTA — 



**E 



to write the reaction cross section as 



(23) 



IV. NUMERICAL RESULTS AND 
COMPARISON WITH DATA 



<Tr=7rA 2 (/o + / 1 +... + / L ) 



7ri?Q— ln<j 1 + cxp 



2.A 2 £ 

1=Il+1,Il+2, 



(l L + l/2)(l L + 3/2)a 



E 



(A0™ +1 d n 
(n + 1)! dl n 



f(l) ■ (24) 



In the special partition by singling out only the lowest 
?i=0 wave for special consideration, then up to the third 
order (AZ) 3 , we obtain 



^In 



l + cxp{e/r} "E 

-— T -Hi) 



1 + exp 



3 a 
4T 



'=1,2,3, 



(25) 



where the derivative in the correction term is 
d 2 



dl 



2 /W 



6(21- 
+ (21 



1 



(?) 


.9 




I (i + s) 2 J 



ON 2 

r 



9 



2 5 2 



(l + .g) 2 (l + 5 )3 



(26) 



Terms on the right-hand side of Eq. (|25|) have direct phys- 
ical meanings. The first term corresponds to the contri- 
bution from the lowest I = partial wave, and the second 
term corresponds to the contribution from I > 1 partial 
waves in the continuum approximation, and the last term 



In presenting our numerical results, we shall label Eq. 
([T]) as Formula I, the sum of the first two terms in Eq. 
(f25"T) [or (|28p ] as Formula II, and the sum of all three 
terms in Eq. (|2"5|) [or (|28p] as Formula III. In simple phys- 
ical terms, Formula I corresponds to the earlier result of 
Ref. Q using the continuum approximation for all partial 
waves. Formula II is obtained by writing out the contri- 
bution from the lowest / = partial wave explicitly and 
treating the higher I > 1 partial wave contributions in 
the continuum approximation. Formula III involves For- 
mula II with the inclusion of corrections up to the third 
order in Al = 1/2. Following Esbensen Q, we shall la- 
bel the cross section obtained in the sum of Eq. over 
the Hill- Wheeler penetration probability, under the as- 
sumption of Eqs. © and (JTUJ), as the Hill- Wheeler cross 
section. 

We examine the sample case for the collision of 
16 0+ 14 N where the experimental data [l3| are shown 
in Fig. 1(a). We show the fit to the fusion cross section 
obtained with E = 8.8 MeV, T = 0.40 MeV (hu = 2.51 
MeV), and ro = 1.45 fm as curves in Fig. 1(a). The differ- 
ences among the three formulas cannot be distinguished 
in the logarithmic plot. In order to see the differences, we 
plot the corresponding ratios of the cross sections rela- 
tive to the Hill- Wheeler cross section in Fig. 1(b). In the 
sub-barrier region, we find that Formula I gives an error 
of about 4.5%, Formula II gives an error of less than 1%, 
and Formula III gives an error of less than 0.01 %. In the 
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FIG. 1. (Color online) (a) The comparison of experimental 
fusion cross section for 16 O + 14 N [llj with theoretical results, 
(b) The ratio of the cross sections from Formulas I, II, and 
III relative to the Hill- Wheeler cross section. 



FIG. 2. (Color online) (a) Comparison of fusion cross section 
data for the collision of identical spin- nuclei 16 + ie O [U, 
HH with theoretical curves, (b) The ratio of cross sections 
from Formulas I, II, and III relative to the Hill- Wheeler cross 
section. 



high energy region, all three formulas give small errors, 
of the order of at most 0.3%. 

We conclude from the results of Fig. 1 that for unequal 
nuclei Formula I is adequate for the sub-barrier region 
if errors of 5% are permitted, Formula II gives a more 
accurate result in all regions with less than 1% error, 
and Formula III gives even smaller errors in all regions. 

In Fig. 2, we examine next the fusion cross section for 
the collision of 16 0+ 16 where the experimental data 
[3 HH are shown in Fig. 2(a) and the theoretical results 
from the three different formulas calculated with the pa- 
rameters E =9.5 MeV, r= 0.4 MeV, and ro=1.3 fm arc 
shown as curves. In this case with identical spin-0 nuclei 
with spatially symmetric wave functions, only the even- 
l partial waves contribute to the reaction cross section. 
Formula II consists of the first two terms on the right- 
hand side of Eq. (|28|) and Formula III consists of all three 
terms in Eq. 

On the logarithmic scale of Fig. 2(a), the results from 
all three formulas cannot be well distinguished. The 
agreement of the experimental data with the theoreti- 
cal curves may appear reasonable. However, in Fig. 2(b) 
we examine the ratio of the cross section obtained with 
the three different formulas relative to the Hill- Wheeler 
cross section. In the sub-barrier region, Formula I gives 
errors of order 6%, Formulas II and III give errors of less 
than 0.3%. 

In the high-energy region, all three formulas give er- 
rors oscillating regularly about zero as a function of E. 
The magnitude of the oscillation is nearly constant for 
Formula II at high energies, but it increases as the en- 
ergy increases for Formula III. These results are in agree- 



ment with the earlier observation of Esbensen [3|], who 
noted that, as a result of the spatial symmetry of the 
wave function such that only even / states contribute, 
the energy separation between the contributing / state 
and the I + 2 state increases as energy increases, and the 
total reaction cross section exhibits a step-wise increase 
when a high-Z state enters into the formation of a fused 
system. As a consequence, the continuum approxima- 
tion contains large and oscillating errors. In mathemat- 
ical terms, the large error arises from the fact that even 
though the expansion parameter Ais = 1/2 is less than 
unity in Eq. (|28[) . it is multiplied by the factor dl v jdv 
with dl v jdv — 2. Thus the effective expansion param- 
eter is (Av)(dli,/dv) = 1 and the expansion in Formula 
III does not properly converge. 

We conclude from Fig. 2 that, for the collision of iden- 
tical spin-0 nuclei at high energies, the continuum ap- 
proximation incurs large errors. Formulas I and II give 
errors of about 5% while Formula III gives even greater 
errors up to 15%. On the other hand, near the sub-barrier 
region Formula II gives very small errors. 

V. BARRIER ANALYSIS FOR 12 C+ l:i C 

The results in the last few sections pertain to the col- 
lisions of both light and heavy nuclei. In the collision of 
heavy nuclei, however, there is the complication that the 
correction term C(l,Ei) for low-l partial waves for the 
barrier analysis must be properly taken into account in 
an iterative procedure, as specified by Eqs. (fT4")) and (fT5"l) . 
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In contrast, for light nuclei collisions, the width parame- 
ter r is found to be substantially smaller than the sepa- 
ration between adjacent barriers so that these correction 
terms can be neglected in the barrier analysis, leading 
to a great simplification of the problem. For simplicity, 
we shall therefore specialize to light nuclei collisions in 
subsequent sections. 

Our ability to reach the simple results in the last sec- 
tions relies on the assumption that the fusion barrier E\ 
for the I partial wave is a linear function of 1(1 + 1), as 
given by Eq. ((9]). There may be nuclear collisions in 
which such an assumption may not be valid. 
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FIG. 3. (Color online) Comparison of fusion cross section 
data for 12 C + 13 C @, [j| with theoretical curves. 
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FIG. 4. (Color online) The dimensionless quantity (oy/VA 2 ) 
as a function of E, on a linear scale (a) , and on a logarithmic 
scale (b) for the collision of 12 C+ 13 C. In the barrier analysis, 
the value of E at which (cr r /7rA 2 ) is 1(1 + 1) + 1/2 is the fusion 
barrier Ei. The horizontal lines indicate the 1(1 + 1) + 1/2 
values shown as numbers given along the right vertical axis. 
The vertical lines indicate the positions of the fusion barriers 
Ei at which a r /irX 2 — 1(1 + 1) + 1/2. Data points are from 
II. 



We examine 12 C+ 13 C where the data are shown in Fig. 
3. The data can be explained well by a coupled-channels 
calculation with the ingoing wave boundary condition 
(IWBC) and the M3Y+repulsive potential ||. Never- 
theless, it is useful to examine these data from a com- 
plementary perspective in barrier-penetrationpoints of 
view. One then finds that the 12 C+ 13 C data d, @ can- 
not be described by Formula I, II or III. Any fit to the 
data near the threshold will miss the data at some other 
energy region. We show in Fig. 3 the results of Formula 
I and II obtained with the parameters Eq = 4.7 MeV, 
r = 0.15 MeV, and r = 1.3 fm as the dashed curve and 
the dash-dotted curve, respectively. The region around 
E ~ 5 MeV is not well reproduced. From the viewpoint 
of barrier penetration and the simple model with the as- 
sumption of Ei = Eq + al(l + 1) in Eq. ([9]), the shape of 
the fusion cross section in the collision of 12 C+ 13 C poses 
a problem. 

To check whether the assumption Eq. (j9]) is valid for 
12 C+ 13 C, we can carry out a "barrier analysis" by plot- 
ting ov/7rA 2 as a function of E as shown in Fig. 4. The 
plots in Fig. 4(a) are on a linear scale and those of Fig. 
4(b) on a logarithm scale. The rule in Eq. (fH]) stip- 
ulates that the barrier Ei is the value of energy E at 
which ov/ttA 2 is 1(1 + 1) + 1/2 + C(l, E{). For light nuclei 
collision for which the width V is substantially smaller 
than the spacing between adjacent barriers, the correc- 



tion term C(l, E{) is small and can be neglected. We plot 
0v/-7rA 2 = 1(1 + 1) + 1/2 as horizontal lines in Fig. 4. The 
energy values E where the horizontal lines meet the data 
points give the locations of the fusion barriers E\ in Fig. 
4. 

On plotting the barrier Ei obtained in the analysis of 
Fig. 4 as a function oi 1(1 + 1), one observes in Fig. 5 that 
Ei is not a linear function of 1(1 + 1), as assumed in Eq. 
©. While the linear relationship is reasonable for I > 2, 
the systematics of Ei for I < 2 appears to be different 
from those with I > 2. 

To examine the problem of barrier penetration, we par- 
tition the partial waves into two parts in the intervals 
[0,Il] and [£j,+l,oo]. We assume the Hill- Wheeler pen- 
etration probability for the partition in [0, Zi], and de- 
scribe the cross section from the partition [Zi+l,oo] by 
the continuum approximation. As given by Eq. (|24[) the 
reaction cross section with the neglect of the correction 
term and the assumption of T; = T is 



2Z + 1 



1=0,1, ...h 



l + exp{(Ei-E)/T} 



ttRI — In <j 1 + exp 



(/ L + l/2)(Z L + 3/2)a 



(29) 



Figure 5 indicates that the dependence of the fusion 
barriers Ei on 1(1 + 1) deviate from a linear relation- 
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FIG. 5. (Color online) The fusion barrier Ei as a function of 
1(1 + 1) in the collision of 12 C+ 13 C. 



ship, for the first three partial waves with I = 0, 1 and 
2. Consequently, we partition the partial waves into two 
partitions of [0,2] and [3,oo] with II = 2. The values of 
Eq,Ei, and E 2 can be read off from Fig. 5 as the starting 
point for parameter search, with minor fine tuning. The 
results with E = 4.7 MeV, Ei = 5.4 MeV, E 2 = 5.85 
MeV, r = 0.15 MeV, and tq = 1.3 fm are shown as the 
solid curve in Fig. 3. We observe that although the fit 
of Eq. ([29]) to the experimental data is not perfect, the 
agreement with experimental data around E ~ 5 MeV is 
substantially improved. The simple comparison indicates 
that a possible solution of the peculiar shape of the fusion 
cross section may involve fusion barriers increasing in a 
non-linear way as a function of 1(1 + 1), corresponding 
to a fusion radial distance occurring at a much reduced 
separation for the lowest partial waves in 12 C+ 13 C col- 
lisions. This may be related to the need for a repulsive 
core in the interaction and the mutual excitation between 
the colliding nuclei as shown in [H . 



VI. BARRIER ANALYSIS FOR THE 
COLLISION OF IDENTICAL LIGHT NUCLEI 

The results in the last section illustrate the applica- 
tion of the rule for barrier analysis in the collision of 
unequal nuclei. The reaction cross section for the colli- 
sion of identical nuclei will need to obey the symmetry of 
the total wave function with respect to the interchange of 
the colliding nuclei. As a consequence, the barrier height 
analysis rule will be modified for the collisions of identical 
nuclei. 

For the collision of identical spin-0 nuclei, the quantity 
(oy/ttA 2 ) is given by 



(oy/71-A 2 ) = 2 



21 + 1 



i=0,2,4,... 



1 + cxp{ (E t -E)] /r,} - 



(30) 



For light nuclei collisions, \Ei — Ei±\\ ^> T/ and we can 
evaluate the above quantity (a r /nX 2 ) at the fusion bar- 
rier Ei analytically. For the collision of identical spin-0 



nuclei, the quantity (cr r /7rA 2 ) at the fusion barrier E\ 
(with even Z value) is 



(oy/71-A 2 ) 



= Z(Z + 1) + 1. 



(31) 



For the collision of identical spin-1/2 nuclei, the total 
spin can be S — or S = 1 , with a weight of 1 /4 and 3/4 
respectively. As a consequence, the reaction cross section 
is given by 

, 1 ^ ~ 1 2Z + 1 



2 x 



- E 

'=0,2,4, 



1 



„ 3 \ - 

2 x - > 

4 ^ 1 

1=1,3,5,7,... 



exp{(£?i - E)\IY X \ 

2Z + 1 
+ cxp{(E l -E)}/T l }' 



(32) 



The barrier analysis rule is different for Ei with even-Z or 
odd-Z partial waves. Assuming \Ei — Ej\ ^> Ti for I ^ j 
for light nuclei collisions, we find that for the collision 
of identical spin-1/2 nuclei, the quantity (a r /n~\ 2 ) at the 
barrier E\ with even Z is given by 



(oy/71-A 2 ) 



1(1 + 1) 



1 

4' 



(33) 



and the quantity a r /-K% at the barrier Ei with odd Z is 
given by 



(ov/ttA 2 ) 



Z(Z + 1) + - 



(34) 



These equations can be utilized to determine the fusion 
barriers from experimental (ov/7rA 2 ) data for the collision 
of identical or non-identical light nuclei. 

We can summarize the rules for the barrier analysis as 
follows. The dimensionless measure of the reaction cross 
section (ay/7rA 2 ) at the fusion barrier E\ is 



K-AA 2 ) 

where K is given by 



= l(l + l) + K, 



(35) 



K = 



^ non — identical nuclei, 

1 identical spin nuclei, 

j even— Z, identical spin | nuclei, 

I odd— Z, identical spin i nuclei. 



(36) 



The differences of (a r /i:X )\e, in the different cases are 
large for the lowest Z = partial wave. The differences of 
(a r /TrX 2 )\E, in the different cases are small, in compari- 
son with the first term 1(1 + 1), when I is large. 



VII. PENETRATION PROBABILITY 
ANALYSIS AND RESONANCES FOR LIGHT 
NUCLEI COLLISIONS 

The penetration probability Pi(E) for the Z partial 
wave is a physical quantity that reveals important in- 
formation on the dynamics of the fusion process. The 
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energy E at which p (E) is 1/2 is at the top of the fusion 
barrier, and the shape of the potential barrier is governed 
by the shape and the energy dependence of Pi{E). It is 
desirable to extract such a quantity from experimental 
data for light nuclei collisions for which the width Y is 
substantially smaller than the separation between adja- 
cent barriers. 



A. Collision of unequal light nuclei 

We shall consider first the collision of unequal light nu- 
clei and express the dimensionless cross section (ov/VA 2 ) 
in terms of the penetration probability Pi (E) as 

= Yl (2* + l)fl(2S) (37) 

™ (=0,1,2,3,.. 

The dimensionless cross section (ov/7rA 2 ) appears so fre- 
quently that it is appropriate to abbreviate it by £(£7) 
that is explicitly a function of the energy E. 

The penetration probability Pi(E) can be extracted 
from the dimensionless cross section £(P) = (ov/VA 2 ) 
if we assume that the contributions of different partial 
waves to the dimensionless cross section arc well sepa- 
rated in energy as in light nuclei collisions. Under such an 
assumption, we can consider the contributions to Pi(E) 
from different partial waves. In the domain of E in which 
Pi{E) is significant, the contribution from each of the 
lower /' < I partial waves is saturated to Pv{E) = 1 
while the contribution form each of the higher V > I par- 
tial waves is negligible. We can decompose the sum over I 
in Eq. (|37p into individual contributions. For the I partial 
wave in the collision of unequal nuclei, the penetration 
probability is then given by 

Pi(E) = E( ^~f (0 9[T(0 - E(25)]e[£(J5)-B(Q], 

(38) 

where is the step function, T(l) is the top delimiter of 
E(-E), and B(l) is the bottom delimiter of £(£')■ For un- 
equal nuclei collisions, the sum of I is over i = 0,l,2,3,.., 
and the delimiters can be shown to be 

T(l) = (l + 1) 2 , (39) 
B(0) = and B(l) = T(l - 1) for I > 1. (40) 

If H(E) = (oy/7rA 2 ) is measured experimentally as a 
function of E, the qpenctration probability Pi(E) for dif- 
ferent partial waves can be determined. 

In Fig. 6 we show the penetration probability Pi (E) as 
a function of E for various I partial waves in the colli- 
sion of 12 C+ 13 C, obtained by using Eq. (|38|) and data 
from [1,0. For a given /, it is possible to determine the 
fusion barrier Ei as the energy at which Pi(E) = 0.5, 
as discussed in an equivalent procedure in Sees II and 
V. We can also extract an empirical width T; where 
2F; is defined as the separation of E between Pi{E) = 



• — • C + C (Notani et al. ('12)) 
*— » 12 C + 13 C (Dayras et al. ('76)) 




2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 

E (MeV) 



FIG. 6. (Color online) The penetration probability Pi(E) as 
a function of E on a linear scale (a) and on a logarithmic 
scale (b) for different partial waves I extracted from the data 
of for the collision of 12 C+ 13 C. 

1/(1 + = 0.731 and P(E) = 1/(1 + e) = 0.269. This 
empirical Yi would be the same as the Yi in the Hill- 
Wheeler formula, if the penetration probability follows 
the Hill- Wheeler formula. 

In Table II, we list the fusion barrier Ei and the width 
Yi extracted from Pi(E) in such a procedure for 12 C+ 
13 C. As one observes, the widths for most of the partial 
waves are about equal to 0.15 MeV except for the I = 1 
and 2 partial waves, which are about 0.11-0.12 MeV. 



TABLE II. The empirical values of E t and Ti from Pi{E), as 
extracted from the data of QjJ for the collision of 12 C+ 13 C. 



1 





1 


2 


3 


4 


5 


Ei (MeV) 


4.91 


5.48 


5.97 


6.43 


6.91 


7.51 


T ; (MeV) 


0.16 


0.11 


0.12 


~0.16 


~0.15 


~0.16 



B. Collision of Identical Spin-0 Nuclei 

We shall consider next the collision of two identical 
spin-0 nuclei. The dimensionless cross section (ov/7rA 2 ) 
in terms of the penetration probability P; (E) is 

^ = 2 m+mi.E). (4i) 

" A / 0.2. U,.. 

Under the assumption that \Ei — Ej\ ^> Yi for I ^ j 
for light nuclei collisions, the contributions of different 
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FIG. 7. (Color online) The penetration probability Pi(E) as 
a function of E on a linear scale (a) and on a logarithmic scale 
(b) extracted from the 12 C+ 12 C data of [H|]. 



partial waves to the dimensionless cross section are well 
separated in energy. For the even-/ partial wave in the 
collision of equal spin-0 light nuclei, the penetration prob- 
ability is given by 

V(E)-B(l) i 



Pi(E) 



2(21 + 1) 



-Q[T(l)-H(E)]e[Z(E)-B(l)], 

(42) 



where we find 



T(l) = l(l + l) + 2l + 2, (43) 
B(0) = 0, and B(l) = T(l - 2) for I > 2. (44) 

The penetration probabilit y Pj (E) extracted from the 
experimental 12 C+ 12 C data and Eq. (|4*2"j) is shown 
on a linear scale in Fig. 7(a), and on a logarithmic scale 
in Fig. 7(b). One finds the fusion barrier for the I = 
partial wave E at 5.46 MeV at which Pi =0 (E) = 0.5. 
The boundary between the I = and I = 2 partial waves 
is approximately at E=5.6 MeV. One observes that Pi (E) 
exhibits resonances. The resonances below E = 5.6 MeV 
are most likely I = resonances whereas those resonances 
above E = 5.6 MeV arc most likely I = 2 resonances. 



C. Collision of identical spin- 1/2 nuclei 

We shall consider next the collision of two identi- 
cal spin- 1/2 nuclei. The dimensionless cross section 
(oy/7tA 2 ) written in terms of the penetration probabil- 
ity Pi{E) is 



E 



{=0,2,4,6, 



3 
'1 



E 



(=1,3,5,7,. 



(2Z + 1)P,(£). (45) 



13 C+ 13 C (Trentalange et al. ('88)) i = 2 



1=4 




2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 
E (MeV) 

FIG. 8. (Color online) The penetration probability Pi(E) as 
a function of E on a linear scale (a) and on a logarithmic scale 
(b) extracted from the 13 C+ 13 C data of [13]. 



Under the assumption that \Ei — Ej\ 3> T; for / ^ j 
for light nuclei collisions, the contributions of different 
partial waves to the dimensionless cross section are well 
separated in energy. In the collision of identical spin-1/2 
nuclei , the penetration probability for the even-i partial 
wave is given by 

P l( E ) = ^l^woW ffl - £(£)]©P(£) - B(l)}, 



(2l + l)/2 



(46) 
(47) 



where T(l) with even I is given by 

T(l) = 1(1 + 1) + (I + l)/2, 

and the function B(l) is given by 

B(0) = and B(l) = T(l - 1) for I > 1. (48) 

In the collision of identical spin-1/2 nuclei, the penetra- 
tion probability for the odd-/ partial wave is 

p i( E ) = W—r^®m) £(£)]ep(£) - B(i)], 



3(2/ + l)/2 



(49) 
(50) 



where T(l) with odd / is 

T(/) = /(/ + l) + (3/ + 3)/2, 
and the function B(l) is given again by 

B(0) = and B(l) = T(l - 1) for / > 1. (51) 
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Using Eq. (gBJl or (gSJ) and the 13 C+ 13 C data from 
[l7| . we extract the penetration probability Pi(E) as a 
function of -E 1 for various I partial waves. The results are 
shown on a linear scale in Fig. 8(a) and on a logarithmic 
scale in Fig. 8(b). In Table III, we list the fusion bar- 
rier Ei and the width Ti extracted from Pi{E) in such 
a procedure for 13 C+ 13 C. As one observes, the widths 
for the even-Z states and the widths from the odd-/ states 
appear to fall into two different groups, with T for the 
even-Z states in the 0.7-0.12 McV range, while the width 
parameters T for the odd-Z states lie in the 0.19-0.22 MeV 
range. There seems to be strong dependence on the even 
or odd property of the angular momentum Z of the fused 
system. 



TABLE III. The empirical values of Ei and Fi from Pi(E), as 
extracted from the data of [TtJ for the collision of 13 C+ 13 C. 



1 





1 


2 


3 


4 


5 


Ei (MeV) 


4.70 


5.59 


6.04 


6.68 


7.16 


7.86 


T ; (MeV) 


0.12 


0.19 


0.07 


0.22 


0.10 


0.20 



It should be noted that in deriving the barrier and pen- 
etration probability rules, we have made the assumption 
that \Ei - Ej\ » T; for Z ^ j. The results in Table II 
and III for £) and indicate that such an assumption is 
substantially valid and is a reasonable and approximate 
idealization for light nuclei collisions . The extracted bar- 
rier height and penetration probabilities are approximate 
quantities that reveal the gross features of the fusion pro- 
cess. 



VIII. CONCLUSIONS AND DISCUSSION 

By treating the angular momentum as a continuous 
variable, the reaction cross section can be evaluated in 
a simple analytical form. The continuum approximation 
of the discrete angular momentum variable carries errors, 
and these errors can be evaluated and amended to pre- 
vious results. 

Three different formulas have been presented in the 
present formulations. Formula I corresponds to the ear- 
lier result of Ref. Q using the continuum approximation 
for all partial waves. Formula II is obtained by writing 
out the contribution from the lowest I = partial wave 
explicitly and treating the higher Z > 1 partial wave con- 
tributions in the continuum approximation. Formula III 
involves Formula II with the inclusion of corrections up 
to the second order in AZ = 1/2. 

For the collision of unequal nuclei, the better formula is 
Formula II, which incurs errors of order 0.7% in the sub- 
barrier regions and errors of order 0.2% at high energies. 
The simpler Formula I incurs errors of about 4.4% in the 
sub-barrier region, and errors of about 0.4% at high ener- 
gies. Higher order corrections in Formula III can be used 



if high accuracy is desired, with errors of about 0.005% 
in the sub-barrier region, and errors of about 0.12% at 
high energies. 

For the collision of identical spin-0 nuclei, the applica- 
tion of these formulas incur substantial errors. The best 
formula for identical spin-0 nuclei is Formula II, which 
incurs errors about 0.2% in the sub-barrier regions and 
errors of about 5.0% at high energies. On the other hand, 
the simpler Formula I incurs errors of about 6.0% in the 
sub-barrier region, and an errors of about 5.5% at high 
energies. 

Simple rules have been presented to determine the bar- 
riers Ei and the penetration probabilities Pi(E) for dif- 
ferent Z partial waves from experimental data, for the 
collision of identical or non-identical light nuclei. The di- 
rect determination of the physical quantities as a function 
of Z gives new insight in the fusion process. The barrier 
analysis rule has been successfully applied to examine the 
relation between the fusion barrier and Z for the patholog- 
ical case of 12 C+ 13 C. The application of the penetration 
probability analysis reveals quantitatively the resonance 
structure in 12 C+ 12 C collisions. 

Wc note that the partitioning of the partial waves into 
the lowest Z region and the higher Z region has some 
advantages in phenomenology. There are situations in 
which the properties of the potential barriers for the low- 
est Z states may deviate from the systematics of those for 
the higher I states. These lowest Z states may need to 
be specially handled. One may provide a different de- 
scription of the penetration probabilities for the lowest 
partial waves, with contributions from higher Z partial 
waves represented analytically in the continuum approx- 
imation. By this partition, the new degrees of freedom, if 
any, can be incorporated into the penetration probability 
to provide a clearer picture of the dynamics of the fusion 
process. 

For simplicity, we have carried out the barrier anal- 
ysis and the penetration probability analysis for light 
nuclei collisions. For collision with heavy nuclei, how- 
ever, r is not small compared to adjacent barrier separa- 
tions \Ei — Ei±i\. The barrier analysis for low-Z partial 
waves needs to be carried out iteratively. While analyti- 
cal expressions have been obtained to carry out such an 
iterative procedure, whether such a barrier analysis for 
heavy-nuclei collisions may be practical remains to be 
investigated. 
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